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We present a model including diffusion and particle-size dispersion for separation of particles 
in deterministic lateral displacement devices also known as bumper arrays. We determine the 
upper critical diameter for diffusion-dominated motion and the lower critical diameter for pure 
convection-induced displacement. Our model explains the systematic deviation, observed for small 
particles in several experiments, from the critical diameter for separation given by simple laminar 
flow considerations. 



I. INTRODUCTION 
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In 2004 lHuang et al.l [l| developed the elegant method 
of particle separation by deterministic lateral displace- 
ment in so-called microfluidic bumper arrays. The 
method, which relies on the laminar flow properties char- 
acteristic of microfluidcs, shows a great potential for fast 
and accurate separation of particles on the micrometer 
scale [H H, H, 0| • Among the key assets of the determin- 
istic lateral displacement separation principle are that 
clogging can be avoided because particles much smaller 
than the feature size of the devices can be separated, that 
the devices are passive, i.e. the particles bump into solid 
obstacles or bumpers, and that the separation process is 
continuous. 

More precisely, particle transport in microfluidic 
bumper arrays is primarily governed by convection due 
to the fluid flow and displacement due to interaction with 
the bumpers in the array [IJ. These processes are deter- 
ministic and the critical diameter for separation of rela- 
tively large particles in these devices is well understood 
in terms of the width of flow lanes bifurcating around the 
bumpers in the periodic arrays However, if bumper 
arrays and particles are scaled down, diffusion will influ- 
ence the separation process and affect the critical parti- 
cle size significantly. Previously reported data on sepa- 
ration of particles in bumper arrays all show a bias to- 
wards larger critical particle size than that given by the 
width of the flow lanes nearest to the bumpers of the ar- 
ray [J, H, Q . In this work we extend existing models by 
adding diffusion and taking particle-diameter dispersion 
into account, and thereby explain the observed discrep- 
ancy. 

In bumper arrays particles are convected by the fluid 
flow through an array of bumpers placed in columns sepa- 
rated by the distance A in the flow direction, see Fig.QJa). 
For a given integer N , The array is made .ZV-periodic in 
the flow direction by displacing the bumpers in a given 
column a distance X/N perpendicular to the flow direc- 
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FIG. f : (a) An array of bumpers (black dots) with the defini- 
tion of the lane number I (small arabic numbers), the column 
number c (large arabic numbers), and the gap number g (ro- 
man numbers). The shaded region illustrates how the shift 
in the position of gap II follows the geometry of the array, 
(b) Close-up of a single gap between two bumpers (disks) in 
the array. Each of the four flow lanes carries the same flow 
rate. Due to the parabolic flow profile in the gap region, the 
width Wj'' of flow-lane I depends on the position in the gap. 



tion with respect to the bumper positions in the previous 
column. Due to this periodicity of the array and the 1am- 
inarity of the flow, the stream can naturally be divided 
into N lanes, each carrying the same amount of fluid flux, 
and each having a specific path through the device, see 
Ref. Q|. 

For a given steady pressure drop, the fluid in the device 
moves with an average velocity uq. Assuming a parabolic 
velocity profile u(x) in the gap of width w g between two 
neighboring bumpers, see Fig.QJb), 



u(x) — 6uq — ( 1 

w s V w„ 



(1) 



the total flow rate Qtot is given by 



?tot = / u(x) dx = w e u . 
Jo 



(2) 



*URL: http : //www .mic . dtu . dk/microf luidics 



For an iV-periodic array, the N flow lanes in a given gap 
carry the same flow rate Qtot/N. The width tu, of lane I 
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is found by solving 



Qt 
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u{x) dx, 



(3) 



where = 



CO 



is the starting position of lane /. 



In the simple bifurcating flow-lane model [J, H[ the criti- 
cal diameter d c is given as d c /2 = w\ . A small particle 
with d < d c will never leave its initial flow lane and will 
thus be convected in the general flow direction follow- 
ing a so-called zigzag path. Large particles with d > d c 
will quickly bump against a bumper and from then on 
be forced by consecutive bumping to follow the skew di- 
rection of the array geometry, the so-called displacement 
path. When a particle gets bumped by a bumper in the 
array it will be displaced perpendicular to the flow di- 
rection until its center is located one particle radius d/2 
from the surface of the bumper. This corresponds to n\ 
lanes of displacement, 
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u(x) dx = N- — - 
Awl 



(4) 



In the bulk fluid, where the lanes are assumed to have 
equal width X/N, see Fig. HJa), the displaced distance 
^disp is therefore 



4isp(rf) = ni— = X 



'Awj 
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(5) 



In this work we extend the simple bifurcating flow- 
lane model by including diffusion and particle-diameter 
dispersion. 
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FIG. 2: O ur model applied to th e exper imental data of 
llnglis et al.l 0| , iHuang et al.l [lj and iLarsenl [f| . Particle di- 
ameter d over gap width w s is plotted as a function of the 
inverse period l/N. The black lines show the theoretically 
pre dicted critic al p article size f or the b umper arrays used 
by llnglis et al] IHuang et al.l [lj and ILarsenl respec- 
tively. The experimental data points are representing par- 
ticles following bumper paths (open symbols), zigzag paths 
(solid black symbols), and neither of these paths (solid gray 
symbols). 



between the bulk displacement and the shift in position of 
the next bumper. The probability p esc for this to happen 
is given in terms of the error function as 



II. MODEL INCLUDING DIFFUSION 

During the average time r = X/u it takes a particle 
to move by convection from one column to the next, it 
also diffuses. We assume that the diffusion process per- 
pendicular to the flow direction is normally distributed 
with mean value zero and variance 



a 1 = 2Dt, 



(6) 



where the diffusivity D is given by the Stokes-Einstein 
expression 



D 



k B T 
Snr/d 



(7) 



Here fee is Boltzmann's constant, T is the temperature 
and r\ is the viscosity of the fluid. 



^disp(^) — Jf X 

V2(j(d) 



(8) 



where we have introduced the <i-dependence explicitly. 
When a particle is transported through a bumper array 
it must bump at every bumper within one period of the 
array in order to be displaced one gap at the outlet. Thus, 
if the particle escapes at least one time in TV attempts, it 
will not be displaced. We define the critical particle size 
d c to be the size for which half of the particles escape 
bumping as they pass one period of the array. Thus d c 
can be found by solving 



N 

E 

k=l 



p esc (d c ) k 1 -p esc (d c ) 



N-k 



(9) 



A. Diffusion model 

In order to escape bumping, a particle must in the time 
interval r diffuse more than the difference €di S p — X/N 



In Fig. we have plotted the critical particle size as a 
function of the bumper period for para meter values c or- 

Hi 

5] 



responding to the bumper arrays used bv lHuang et al 
(dashed line) . llnglis et all @ (dotted line), and Larsen 
(full line). 
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B. Comparison to experiments 

The observation that the critical particle size in a 
bumper device is larger than the width of the first flow 
lane is also supported by the experimental data in Fig. 2 
of Ref . [3j . Our model suggests that the deviation of the 
critical particle size from the width of the first flow lane 
can be explained by diffusion of the particles. In Fig. [2] 
above, it is seen how well the lines are predicting the 
transition between zigzag paths and displacement paths: 
the thick line is dividing open and closed triangles, the 
dashed line is dividing open and closed squares, and, to a 
lesser extent, the dotted line is dividing open and closed 
circles. 

Using parameter values corresp onding to the bumper 
device presented by iHuang et al.1 [l[ our model predicts 
a critical particle diameter of 0.45 times the width of the 
gap for the particles traveling through the device with an 
average velocity of 400/im/s. This is in good agreement 
with Fig. 2(a) in Ref. fl. 



III. A DISCRETE MODEL INCLUDING 
DIFFUSION AND DISPERSION 

The particles used in the experiments on separation of 
particles are not mono-disperse. Their average diameters 
are distributed around a certain mean value with a rel- 
ative standard deviation Ad/d, which typically is 20 %, 
10% and 5% for particles with d — 25 nm, d = 100 nm 
and d = 500 nm, respectively. 

In the following we introduce a discrete model of the 
transport of particles with different diameters d in an N- 
periodic bumper array taking convection, diffusion and 
sizc-dispersion into account. The model allows us to 
study the relative influence of all three phenomena on 
the separation efficiency in a fast and simple manner. 
In particular our results suggest that the critical size for 
separation or displacement, studied above, must be sup- 
plemented by a smaller critical size below which pure dif- 
fusion governs the motion of the particles in the bumper 
array. This prediction has not yet been tested experi- 
mentally. 



A. Definition of the discrete model 

At any instant, a particle is assumed to be positioned 
in the center of a flow lane I situated in a gap g in some 
column c of the array. For simplicity we further assume 
that the size distribution of any given set of particles is 
a normal distribution with a mean value given by the 
size quoted by the manufacturer and a relative standard 
deviation of 10 %. 

By convection any given particle moves from one col- 
umn to the next. If it ends up in the last lane (I = N — 1) 
in one gap, it will be shifted to the first lane (I = 0) in 
the subsequent gap. Otherwise it will stay in the current 



gap and move up one lane. In our model pure convection 
is therefore described by the discrete map 

(c,0,i) | -M, , , 1N , . ( 10 ) 
I (c + 1, g, I + 1) , otherwise. 

Because of the finite diameter d of the particle there is 
a minimum and a maximum lane number that it can oc- 
cupy. The minimum lane number is the smallest integer 



w{ > -. 

1=0 



(11a) 



Similarly, the maximum lane number / max is the largest 
integer that satisfies 
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(lib) 



Consequently, the simple convection mapping from 
Eq. (fTOl needs to be modified to account for the finite 
size of the particles 

(c + l,g+ 1,U„) ,if/ = iV-l, 
(r. lh l) - { (c+l,g,l + l) , if I < Z max - 1, (12) 
(c+1 , g, ^max) , otherwise. 

The above convection scheme accounts for the separation 
of particles in deterministic lateral displacement devices 
according to size. The critical particle radius predicted 
by this model is 



2w 



(0) 



(13) 



in accordance with the geometric arguments of Ref. Q . 

To characterize the separation quantitatively, we define 
the relative particle numbers rg, r\ and as 

r = relative number of particles following the (14a) 
zigzag path. 

r\ = relative number of particles following the (14b) 
displacement path. 

rd = relative number of all other particles. (14c) 

With these definitions the sum ro+ri+rd is always unity. 
If t*o = 1 all particles follow the zigzag path and if r± = 1 
all particles follow the displacement path. If rd 7^ some 
of the particles end up at positions not explained by the 
deterministic analysis of the separation process. In Fig. [3] 
we have plotted the relative particle numbers ro, r\, and 
rd as a function of average particle diameter d. 



B. Pure mono-disperse convection 

For mono-disperse and non-diffusing particles, the 
model results, as expected, in two modes: the zigzag 
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FIG. 3: Reiative number ro, ri and ra of particles following the zigzag path, the displacement path, and neither of these two 
paths, respectively, plotted as a function of the normalized, average particle diameter d/d c . The parameters of the bumper 
array is taken form Ref. [f|: N = 100, A = 8 /urn, w s = 1 (im, L = 207V A = 16 mm, and uo = 250 pim/s. The buffer liquid is 
water at room temperature. Neglecting diffusion (solid symbols), the particles follow the zigzag path if d < d c — 118 nm. If 
d > d c they follow the displacement path. Including diffusion (open symbols), the small particles are not influenced by the 
bumpers. For d > 2.1 d c the influence of the bumpers set in, and for d > 3.4 d c the particles follow the displacement path. 
The thick black vertical lines indicate the particle size when small particles stop behaving purely diffusive (left-most lines) and 
when large particles begin a purely deterministic displacement (right-most lines). 



mode and the displacement mode, see the closed circles 
in Fig. [3] For d < d c we have ro = 1, and for d > d c we 
have r\ = 1, while we always have = 0. The relative 
particle numbers can therefore be written as 



(ro,r d ,n) 



(1,0,0) ford<d c 
(0, 0, 1) for d > d c 



(15) 



C. Influence of size dispersion 

If we assume that the particles are not mono-disperse, 
but are distributed around a mean size d with standard 
deviation Ad, the shift as a function of d from the zigzag 
mode to the displacement mode happens gradually in- 
stead of abruptly at a certain critical size d c (Fig. [3l 
closed squares) . The relative number of particles follow- 
ing the zigzag path ro can be found by integrating over 
all particle sizes smaller than the critical diameter given 



by the array geometry 
1 



y/2TT(Ady 



:exp '^F> ds - (16a) 



Similarly, the relative number of particles following the 
displacement path r\ can be found by integrating over 
all particle sizes larger than d c 



ri 



1 ( ~{s-df 
^(Ad) 2 6XP V 2(Ad)2 



ds. (16b) 



The system is still a bi-modal system because = for 
all particle sizes. 



D. Influence of diffusion 

In ID during the time step r a particle diffuses the 
distance i, the average of which is the size-dependent 
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diffusion length <j(d) given by 



a(d) = {£) = \Z2Di 



, k B T A 

37T77<i Mq 



(17) 



In our model we discretize the transverse diffusion as the 
properly rounded number of flow lanes crossed by the 
particle during diffusion, 
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(18) 



The addition of diffusion smears out the displacement of 
the particles and causes the critical radius to be larger 
than in the diffusion- less case (Fig.[3j open symbols). 



1. Bumping criterion for small particles 

Very small particles are completely dominated by dif- 
fusion, and the particle distribution at the end of the 
array is simply given by the transverse diffusion of the 
particles during the time L/uq it takes for the particle to 
be convected all the way L through the array. For small 
particles we therefore have 
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(19) 



where a 2 = 2DL/uq. 

As the particle diameter d is increased, the bumpers 
begin to become important as the diffusion length (r(d), 
Eq. (fT7|) , is decreased and becomes equal to the displace- 
ment length £dis P , Eq. ©. Using the criterion a(di) = 
^disp, with the parameter values used in Fig. [3l we find 
that particles stop behaving as small diffusion-dominated 
particles and start interacting with the bumpers when 
d\ = 2.1 d c . This cross-over value is indicated by the 
left-most vertical lines in Fig. [3l and it fits well with the 
simulation data. 



2. Bumping criterion for large particles 

Large particles will interact with the bumpers at every 
row in the array and their position is thus reset at every 
bump to be ^disp- Diffusion can therefore be neglected 
for such particles, they all follow the displacement path, 
and n = 1. 

As the particle diameter is lowered, the probability p e sc 
that a particle escapes the displacement path can be esti- 
mated by the probability of diffusing from the displaced 
position to the last flow lane in the gap, i.e. the distance 
^disp - jj- This probability p csc is given by Eq. ©. 

In order to follow the displacement path, a particle 
must bump at each row in the array. If we consider a N- 
periodic array with m full periods, the particles will have 



mN bumping opportunities as they pass through the en- 
tire array. If a particle evades bumping at a bumper, it 
will be convected by the flow through a full period of the 
array before bumping is possible again. Because of the 
escape, it will miss N bumping opportunities and end 
up one gap from the displacement path. Consequently, 
if a particle escapes one time, it will only have (m — 1)N 
bumping opportunities and has therefore escaped bump- 
ing with a probability of l/[(m — 1)N]. The upper crit- 
ical particle size di for convection-induced displacement 
is defined using Eq. §8§ as 



p esc (d 2 ) 



1 



(m - 1)N 



(20) 



For the device used in the experiments by ILarsenl [f| we 
find g?2 = 3.4g? c . The predicted upper limit for diffusion 
dominated motion and the lower limit for convection- 
induced displaceme nt comp ares well with the experimen- 
tal observation by ILarsenl [j| that some particles end 
up between the displacement path and the zigzag path 
(Fig. [2l gray triangles) . 



IV. CONCLUSION 

Experimental data on separation of particles in 
bumper arrays all show a systematic deviation from pre- 
dictions made from the bifurcating flow lane model [3, 
d, [H. By adding diffusion to the bifurcating flow lane 
model, we explain the discrepancy. In addition, we have 
suggested a simple discrete model for simulating particle 
separation in bumper arrays. The presented model takes 
particle diffusion and size dispersion into account and has 
been validated against experimental data for a bumper 
device with period N = 100. The transition from zigzag 
path to displacement path happens at a particle size 2.1 
to 3.4 times larger than the critical particle size predicted 
from geometrical arguments. This is in correspondence 
with the experimental observations from Larsen Q . Our 
discrete model and the estimates presented in this paper 
suggest that particles of twice the size of the geometrical 
critical size of the N = 100 bumper device behave diffu- 
sively and are not affected by the bumpers because the 
small diffusive particles rarely come into contact with the 
bumpers due to random Browninan motion. 
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